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Top-antitop pairs produced in the decay of a new heavy resonance will exhibit spin cor- 
relations that contain valuable coupling information. When the tops decay, these correla- 
j-Gj" tions imprint themselves on the angular patterns of the final quarks and leptons. While 

many approaches to the measurement of top spin correlations are known, the most com- 
J> , mon ones require detailed kinematic reconstructions and are insensitive to some important 



CO 



spin interference effects. In particular, spin-1 resonances with mostly- vector or mostly- axial 



couplings to top cannot be easily discriminated from one another without appealing to 
^j' ■ mass-suppressed effects or to more model-dependent interference with continuum Standard 
Model production. Here, we propose to probe the structure of a resonance's couplings to 
tops by measuring the azimuthal angles of the tops' decay products about the production 
axis. These angles exhibit modulations which are typically 0(0.1-1), and which by them- 
^ \ selves allow for discrimination of spin-0 from higher spins, measurement of the CP-phase for 
spin-0, and measurement of the vector/axial composition for spins 1 and 2. For relativistic 
tops, the azimuthal decay angles can be well-approximated without detailed knowledge of 
the tops' velocities, and appear to be robust against imperfect energy measurements and 
neutrino reconstructions. We illustrate this point in the highly challenging dileptonic decay 
mode, which also exhibits the largest modulations. We comment on the relevance of these 
observables for testing axigluon-like models that explain the top quark Afb anomaly at the 
Tevatron, through direct production at the LHC. 



I. INTRODUCTION 



New resonances that decay to it are a common feature of many models of physics beyond 
the Standard Model (SM). Indeed, top's unusually large coupling to the electroweak sector 
is often taken as a hint of some special role in TeV-scale physics. 

Recently, another hint has arisen at the Tevatron, in the measurement of a 3.4a anomaly 
in the top pair forward-backward asymmetry (Aps) at high invariant mass [l| (as well as 
lower-significance biases of the same sign observed in and more recently in the dilep- 

tonic mode fl). One possible explanation of this effect is the presence of an s-channel 
qq — > it resonance somewhere beyond the Tevatron's direct discovery reach, which is inter- 
fering with the SM production. As argued in [6], a straightforward model framework which 
can lead to the required interference is a color octet spin-1 boson with sizable axial couplings 
(of opposite sign) to light quarks and tops; e.g., some variant of a heavy axigluon [7H9|. 

Such a resonance should be visible at the LHC, even if it turns out to be very broad, and 
likely even in the upcoming 7 TeV run This discovery would already be extremely 

suggestive. Direct confirmation of Ap^-related effects may also be pos sible, despite the 
fact that the pp initial state is forward-backward symmetric ll|, [l4- 17]. However, fully 
establishing the connection between the resonance and Afb will benefit from independent 
measurements of the resonance's chiral couplings to top quarks and light quarks, which 
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will necessarily involve additional detailed kinematic analyses both on- and off-peak 
Most of these analyses exploit effects that, like Aps itself, arise from interference with QCD. 

Motivated by this endeavor, but also with an eye toward more general models, we can 
consider the following question: Given that a resonance (either a peak or high-mass excess) 
is discovered in the m t i spectrum at the LHC, how much can we directly learn about the 
chiral structure of its couplings to tops without relying on more model-dependent QCD 
interference effects? In particular, in the case of a spin-1 resonance, could we determine if 
it couples to tops mostly-axially by using only on-peak observables? 

Of course the presence of a large axial component could be readily established if the 
individual chiral couplings are highly biased toward left-handed or right-handed. At energies 
well above 2m f , the different chiralities map almost uniquely onto distinct top helicities, 
which can be discriminated by studying the rest frame decay angles of individual tops 231.124 1. 



The differences are also encoded in the lab-frame energies of the leptons and fe-jets relative 
to the top and to each other 25M271]. as well as in the energy sharing between the jets (or 



subjets) in hadronic decay 28]. 

The situation becomes more subtle if the couplings to the two top chiralities are more 
democratic. Single-side spin analysis then becomes inadequate to accurately probe the rela- 
tive magnitudes between vector and axial couplings, and we must incorporate measurements 
of spin correlations between the two sides of the event. These correlations are affected by 
the interference between the different helicity production modes, and can be used to directly 
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probe the signed ratio of top's chiral couplings. Nonetheless, the most commonly used top 
spin-correlation variables 29-34J (reviewed in 23, 35 1) are largely insensitive to helicity in- 
terference for heavy spin-1 resonances, and generally involve complicated reconstructions 
that are not obviously ideal for fast-moving tops. Our main aim in this paper is to demon- 
strate how to construct simple observables that are sensitive, as well as straightforward to 
measure. 

For spin-1 resonances well above 2mt, the on-peak helicity interference is most directly 
encoded in the azimuthal distributions of the tops' decay products about the it production 
axis. This is quite analogous to the effect measured at LEP in Z — > r + r~ |36l-l39j] . It is also 
closely related to the observables explored in 4CH44| , which use the same type of interference 
as a probe of the spins of new pair-produced particles, or particles that are singly-produced 
in association with a jet. We now assume that we know the spin of the top quark, and use the 
interference to tell us about how it was produced^ As is typical for top spin observables, the 
decay product most sensitive to the correlation is the charged lepton, and we focus here on 
the dileptonic decay mode. We also emphasize that the same effects occur in the £+jets mode, 
though with smaller amplitude. We will show that the dileptonic azimuthal distributions 
can be measured at the LHC, even given the usual kinematic difficulties, and even using 
very minimalistic reconstruction strategies. These distributions exhibit surprisingly robust 
modulations at the 30-40% level for pure vector or axial coupling, but with opposite signs. 

Our observations also generalize beyond spin-1 resonances. The minimal couplings to 
spin-2 have a similar chiral structure, which will be amenable to the same style of analysis. 
Spin-0 resonances near it threshold have been a topic of detailed study for the past two 
decades, due to the possibility of detecting scalars from the Higgs-sector via gluon fusion 



[e.g. [46|, |47j). Variables sensitive to the scalar/pseudo-scalar composition have therefore 
been known for a long time. However, we will emphasize that azimuthal correlations make 
this composition visible in a very transparent way, especially if the resonance is heavy. They 
also allow for a straightforward determination of whether a resonance is spin-0 versus higher 
spin, independent of it polar decay angle correlations or the production angle distribution. 
Given the large amount of easily-measured information made available to us through these 
correlations, we highly advocate incorporating them into the set of it resonance discrimina- 
tion observables laid out in 



23]. 



In the next section, we discuss the angular spin-correlation observables. We start by 
reviewing the basic formalism, and then proceed to the details particular to spin-0, spin-1, 
and spin-2. In particular, our discussions of spin-1 and spin-2 demonstrate how vector- 
coupled and axial-coupled resonances can be distinguished from each other with azimuthal 
correlations, as well as from spin-0. In Section IIHI we demonstrate that our observables 
should be^ccessible experimentally in the dileptonic channel. We estimate the amount 

for a study using the individual top decay azimuthal angles to probe new it production 



Sec also 
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contributions at an ILC. 
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of cross section required to be able to distinguish between pure vector and axial-vector at 
early- and late-stage LHC, and discuss possible relevance to the Tevatron top Afb anomaly. 
We conclude in Section IIV1 The Appendices include complete formulas for leading-order, 
on-peak ti production and decay angular distributions for resonances up to spin-2. 



II. SPIN CORRELATIONS 

After a heavy resonance decays to a pair of top quarks, the tops themselves decay into 
b£ + v or bqq' through an on-shell W boson. In general, a given 6-body final-state kinematic 
configuration could be produced by a number of intermediate top/anti-top spin configu- 
rations. Since we never measure the tops' spins, these different spin channels coherently 
interfere with one another, leading to distinctive imprints on the final-state angular distri- 
butions. 

The relevant behavior of the top decay matrix elements becomes transparent when we 
exploit a fortuitous simplification. If we look at the 3-body semileptonic top decay amplitude 
at leading-order, stripping off the W propagator, we can Fierz down to 

M(t^b£ + u) ex (u{t) T L ev{£) L ) (u(6)£ eu(z/)l), (1) 

where w's and i>'s represent the spinor polarizations, and we have kept only the two compo- 
nents from left-chirality. Given a specific top spin and decay configuration, we can map out 
how the amplitude changes as we rotate the decay products in the top rest frame. Because 
the 6-quark and neutrino factorize off into a rotationally-invariant product, the angular 
dependence of the amplitude is determined exclusively by the direction of the lepton: 

M(t t ^b£ + u) oce^ /2 cos^ 

M(t^ -> btv) oc e-^ /2 sin |, (2) 

where we are using a standard right-handed coordinate system with spin quantized along 
the z-axis. If we instead decay an anti-top, we simply swap cos <H> sin. (The same formulas 
also hold with £ — > d/s in fully hadronic decay.) Notably, the phases of the amplitudes 
depend on the azimuthal angle of the lepton about the top's spin vector. 

Since we are ultimately interested in resonances beyond the Tevatron reach, and well- 
above top threshold, we will restrict our discussions to the chiral production limit. (More 
complete formulas, valid down to threshold, appear in the Appendices.) In this limit, the 
most natural spin basis is the helicity basis. We therefore quantize spin along the top/anti- 
top production axis in the ti CM frame, and measure around this axis (the "azilicity angle," 



in the terminology of 48]). By convention we choose +z to point along the t direction. This 
may differ with conventions used in other places, where the i decay coordinate system is 
independently defined with respect to its own direction of flight. To define the x- and y-axes, 



3 



we use the beams. We can construct y by crossing the t direction with one of the beam 
directions, so that it points out of the ti production plane. (The ambiguity over which beam 
to pick at the LHC will not be important for our purposes.) The x direction is then defined 
as usual for a right-handed coordinate system, as y x z. Consequently, when we measure <ft 
angles about the production axis, we define <fi = to lie in the production plane. This is, of 
course, our only physically-motivated choicej§ Finally, we construct the individual top and 
anti-top decay coordinate systems by actively boosting the tops to rest (without rotation) 
along the production axis, and measuring all angles in this common reference frame. 

Assembling a complete 6-body decay matrix element for a hypothetical heavy particle 
("X") is now straightforward. Taking the narrow width approximation, and implicitly 
picking some definite initial spin state, we can write 

Mot = J2 M ( X ~+ tata)M(t a ~> b£ + u)M(i a -> b£~u), (3) 
a. a 

where a and a label the spin of the top and anti-top, respectively. For spin-0 resonances, 
the tops are forced to have opposite spins (identical helicities). In the chiral limit, tops from 
a heavy spin- 1(2) decay will have identical spins (opposite helicities). The matrix elements 
schematically reduce to 

M tot (X -»• btvbi-v) ~ M(X -> t t tj e'fo-t'V 2 cos ^ cos ^ + 

M(X -+ t&) e -*te-^)/ 2 sin | sin | 
M tot (X m -+ btvlrv) ~ M(X 1{2) -+ t t t t ) jfa+W cos | sin | + 

M{X l{2) -+ t&) e<**#M* sin | cos |, (4) 

with the bar indicating the anti-top (£~) decay angle. The final rate will therefore contain 
modulating terms dependent on (^±<^, with phase or magnitude dependent on the complex 
production matrix elements. For the case of spin-0, which modulates as (fa — (f> e ), the two 
X decay elements are conjugates of each other. Consequently, the modulation is subject to 
a pure phase offset, where phases differing from or ir signal CP- violation. For spin-1 and 
spin-2, which modulate as (fa + fa), the effect is instead pure magnitude, varying with the 
signed ratio of chiral couplings. (E.g., in the case of purely LH or RH chiral production, 
there is only one helicity channel, and therefore no modulation effect.) Since we will consider 
scenarios where X is a singly-produced resonance, and not necessarily extremely narrow, 
there may also be contributions to the azimuthal modulations from interference with the 

2 In the presence of nonvanishing ti system pj>, we must be somewhat careful when performing this coor- 
dinate construction, as the beams are not necessarily aligned in the it CM frame. We address this issue 
m Section llnl 
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SM. Still, the effect tends to be subleading for resonances with a well-defined peak. We 
discuss the issue in more detail at the end of Appendix [B] 

We can also see similar correlations between any two decay products, one from the top 
side and one from the anti-top side. While Eq. [T] naively suggests no dependence on the b- 
quark and neutrino directions, these are still highly kinematically correlated with the lepton 
via the (approximately) fixed mass of the W boson and overall four-momentum conservation 
for the top. Integrating out the phase space of all but one particle on each side of the event 
leaves over angular correlations that are essentially identical in structure to the double-lepton 
correlations, including phase information. The correlations simply scale with constant spin 
analyzing powers associated to each particle species: 1 for leptons, -0.3 for neutrinos, -0.4 
(+0.4) for fe-quarks (IV-bosons), and 0.5 for the softer of the two W decay products in the 
top rest frame 49[ . The effects which we seek to exploit are maximal for leptons, but could 
also be measured in principle by including fe-jets or jets from the W decay. (Indeed, formulas 
analogous to Eq. HI but instead using H^-bosons, appeared long ago in [31[ in the context of 
spin correlations between top pairs at the SSC.) 

Here, we focus on the double-lepton correlations, and therefore the dileptonic tt decay 
channel. This maximizes the azimuthal modulation effects and is the least susceptible to 
reducible backgrounds. However, it also has the lowest rate (5%, versus 30% for £+jets), and 
suffers from kinematic reconstruction difficulties due to the two approximately back-to-back 
neutrinos. We address the latter difficulty in Section IHIl Incorporating the £+jets mode 
will also be very important, even though the modulations will be smaller in amplitude (e.g., 
40% if we correlate a lepton with a 6-quark or jy-boson from the other side of the event). 
We reserve a detailed exploration for future work, but include some commentary. 

In what follows, we will usually be integrating out all but one or two angular variables, 
only restoring the full (and usually very complicated) correlated angular dependence in the 
Appendices. Much of the discussion will parallel that of [23|, which comprehensively studied 
some of the standard correlation observables in the context of resonance discrimination 
for arbitrary spins. We also note that more sophisticated analyses, such as the likelihood 



approach of 50j, |5jj, could potentially distinguish the different resonance scenarios even 



more efficiently. 



A. Spin-0 

We begin with the well-studied example of a spin-0 resonance. Single-production of 
such a resonance is typically chirality-suppressed from a qq initial state, and dominantly 
proceeds through heavy quark loops in gluon fusion. It is therefore perhaps nongeneric for 
this process to have a healthy rate above backgrounds when the resonance mass is at or 
above a TeV. However, independent of any model-building concerns, at a minimum, spin-0 
serves as both an illustrative example of our methods and as a baseline hypothesis against 



5 



which any higher-spin resonances must be tested. 

At dimension- four, an arbitrary scalar coupling to tt is made of a linear combination of 
CP even and odd pieces, 

4* = -^? {y + iyj 5 )T%V, (5) 

where Tg- is the color matrix appropriate to the scalar's representation (singlet or adjoint). 
The dependence on the scalar's color factorizes, so we ignore it in the remainder of the 
discussion. Since our interest is in the different top helicity components, we can rewrite 
Eq. in the chiral field basis, 

Ant ^-V(f> (e ia i R t L + e- ia i L t R ) , (6) 

where ye ia = y — iy. Since the ti pair can have no angular momentum along their pro- 
duction axis, only same- helicity tops are created in the scalar decay. In the chiral limit, 
the production matrix elements differ only by the overall phase factor e ±ia . The angle a is 
equal to for a CP-even scalar and 7r/2 for a pseudoscalar. 

Squaring .Mtot in Eq. [3] and integrating over phase space except for lepton orientations, 

oc 1 + cos 9g cos 9g — sin 9g sin 9g cos — 4>e + 2a) . (7) 



(Kip dVt 



Getting the distribution of polar angles is trivial from this equation, as either lepton's 
azimuthal integral will kill the last term. This leaves over the well-known formula 

d 2 T 

— =- oc 1 + cos9 e cos9 e . (8) 

a cos 9 1 a cos 9i 

The prefactor of the double-cosine term is characteristic of a scalar, and has been advocated 



for distinguishing the overall spin of a resonance |23| (though, we re-emphasize, typically 
with a coordinate change 9i —> n — 9g). However, we see that dependence on the phase of the 
coupling has dropped out. To retain it, a traditional approach is to instead integrate out all 
angles except for the 3D opening angle between the leptons, which we call Recall that we 
have actively boosted both tops along the production axis into their respective rest frames, 
and then defined a common coordinate system with the z-axis given by the t direction. It 
is in these coordinates where we measure x- We get a simple linear dependence on cos x, 

dT 1 - 2cos(2a) 

ocl + ^^cosx- (9) 



dcosx 



A CP-even scalar has a different slope compared to a pseudoscalar, — | vs. +1, making 
the distinction between the two straightforward in principle, provided that one can fully 
reconstruct the individual top quark rest frames. Such a reconstruction can be challenging. 



Elsewhere in the literature this has been called 6. We relabel it to avoid confusion. 
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However, as shown in 52j, the resolution on the opening angle may still be good even for 
highly boosted dileptonic top pairs. A more fundamental issue, as we will see in the next 
subsection, is that spin-1 and spin-2 resonances all look identical to the CP-even scalar (and 
to each other) when using this angle. 

A more straightforward way to maintain the a-dependence, without needing polar angles, 
is to integrate Eq. [7] over all angles except for the offset between the leptons. This gives 
a modulating distribution 

dT /7T\ 2 

ex 1 - - cos(4> e -4> e + 2a), (10) 



d((j)e-(f)e) V4 

with a sizable amplitude of roughly 60%. We see that the modulation flips sign as we go from 
CP-even scalar (a = 0) to pseudoscalar (a = vr/2), and that the magnitude of the modulation 
is the same 0(1) size for any admixture. Since <p£ and <f)£ are boost-invariant along the 
individual tops' directions of motion in the it CM frame, using Eq. [10] may be amenable to 
simpler ti reconstruction methods that do not necessarily need accurate information about 
the tops' velocities. Moreover, as we will see below, the equivalent distributions for spin-1 
and spin-2 resonances are flat, breaking the degeneracy with the CP-even scalar, and already 
allowing a significant degree of spin discrimination in this variable alone. 



B. Spin-1 

We can write out the generic dimension- four coupling between a spin-1 resonance and a 
ti pair as 

£ int = A^f (g v + 9A1 5 ) (11) 

where gy and g^ are real numbers. Again, we will suppress color in what follows, and go 
over to the chiral field basis, 

Ant A„ (g L i Ll H L + g R i R -ft R ) ■ (12) 

Unlike the couplings of a scalar, these couplings preserve chirality, and they differ in magni- 
tude and sign instead of being complex conjugates of each other. For convenience, we define 
a coupling angle £, 

g L = g cos f 

g R = gsin£. (13) 

Our decay angular distributions will be functions of £, with the overall g factored out. By 
imposing symmetry under parity, we can have the special cases of = g R (true vector, 
£ = 7r/4) or gL = —g R (axial- vector, £ = 37r/4). We are particularly interested in finding 
observables that discriminate between these two cases. 
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We again assume the chiral limit for top production, setting m t = in the decay matrix 
element of the heavy resonance. The finite-mass effects in the angular distributions can 
in fact be different for the true vector and the axial-vector, in principle providing addi- 
tional means of discrimination if the resonance is light enough^) However, by focusing on 
the features of the angular correlations that are asymptotically insensitive to m t /M, our 
observations here will be robust up to arbitrarily high resonance masses. 

We begin our discussion of spin-1 angular distributions by looking at the standard double- 
cosine distribution, integrating over all other angles, 

d 2 T 

=- oc 1 + cos(2£) (cos 6^ — cos 6^) — cos6>£ cos 6^. (14) 

a cos 0£ a cos 0£ 

This is the direct analog of Eq. [S]in the spin-0 case. The first term signals parity violation, 
and is only active if the resonance is neither pure vector nor pure axial-vector. (Note that 
cos(2£) = (g\ — g R ) / (g\ + gj^).) As discussed in the introduction, a large chirality bias could 
be established already with single-side angular distributions, such as dT /d cos 9e, and the 
construction of the double-cosine distribution then essentially serves only as a cross check. 
Incidentally, this observation by itself would already disfavor a spin-0 interpretation. On the 
other hand, in the more chirality-symmetric case (cos(2£) ~ 0, \gi\ ~ \9r\), the correlated 
distribution becomes more useful, as the sign of the double-cosine term is opposite that 
for the scalar. This is a direct indication that the tops are coming out in a same-spin 



(opposite- helicity) configuration, again disfavoring spin-0 |23|. Crucially, however, all terms 
in the distribution are insensitive to the relative sign between the chiral couplings, up to 
the finite-mass corrections mentioned above. Consequently, it will likely be quite difficult to 
tell whether the resonance is vector-coupled or axially-coupled using only the double-cosine 
distribution. 

Turning to the distribution over the 3D lepton angle x does not help. It is 

oc 1 cosx, (15) 



d cos x 3 

independent of £. This is identical to the CP-even scalar distribution, Eq. HI again up to 
the small (mt/M) 2 corrections which we have neglected. We therefore learn very little from 
this distribution if we are interested in the properties of spin-1 resonances. 

To robustly determine the signed ratio of chiral couplings, we again suggest the alternative 
tactic of focusing on the lepton azimuthal angles with respect to the production axis. There 
are now three major differences with the spin-0 case. First, as indicated in Eq. HJ the 

4 In particular, the double-cosine distribution for spin-1 resonances in the it threshold limit is quite different 

between the two cases. An axial- vector has a relative double-cosine coefficient of -1 for all values of m t /M, 

whereas a true vector's relative coefficient approaches -1/3 at threshold (c/. Eq. IB3|) . In contrast, for a 

1 TeV resonance, the coefficients differ from each other by only about 10%, and the discrepency shrinks 

like (mt/M) 2 . Similar conclusions also hold for the cosx distribution, though the finite-mass effects there 

might be easier to measure since the distribution is one-dimensional instead of two-dimensional. 
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modulation switches from (<fig — <pi) to (<fre + 0^), and therefore now refers to the orientation 
of the production plane (zx-plane)Jfl This is a consequence of the fact that the spin-1 
resonance carries a polarization which is necessarily fully correlated with the beams. The 
second difference is that the coupling angle (a — > £) now affects the modulation's amplitude, 
not its phase. The third difference, which also follows from nonzero polarization, is that the 
decay angle distributions can become correlated with the ti production angle, which we 
call G. In what follows, we will implicitly average between rates at O and ti — O, so that 
effectively 6 G [0,tt/2]. 

For a resonance produced in qq annihilation, the component of the resonance's spin along 
the beam is Jbeam = ±1, and the joint production/azimuthal angle distribution is 

m^ITT 1 n K i 1 + cos2 6 ) - (t) 2 sin ^) sin2 9 cos ^ + ( 16 ) 
d((pi + (pi) a cos 6 V4/ 

Unlike the double-cosine distribution, which depends only on cos(2£), the azimuthal mod- 
ulation depends on sin(2£) = Ig^g^jiej^ + g 2 ^) = 2 /(gh/gn + 9rI9l)-, and hence is directly 
sensitive to the signed ratio g^/ 'g^. Perhaps fortunately, the modulation is largest for cen- 
tral production, which is the easiest case to reconstruct in real detectors. (The modulation 
shuts off for ti produced along the beam, since only one helicity state can contribute.) Fully 
integrating out the production angle, the distribution reduces to 



dT Jh =±1 _ 1 /7T\2 

,J beam - J -- L 

d((j>t + 4>i) 



oc 1 - - (^) 2 sin(20cos(^ + ^), (17) 



which is very similar to the spin-0 <p- difference distribution, Eq.[10l However, it is straightfor- 
ward to see that, respectively, a spin-1 resonance is non-modulating over (0^ — 0^)E whereas 
a spin-0 resonance is non- modulating over (0^ + 0^). Using both of these distributions can 
therefore tell us quite a lot about both the spin of the resonance and its coupling structure, 
independent of any other measurements such as the production angle distribution (trivial for 
spin-0 but nontrivial for spin-1), or observables incorporating the tops' polar decay angles 
(sensitive to chirality biases). 

We can also give the analogue of Eq. [16] for a resonance with Jbcam = 0. Such a state 
would come from gluon fusion and is typically suppressed, partly because the gluon parton 
luminosities drop very quickly with energy, and partly because the required gauge-invariant 
operators are dimension-sixjj We include this case for completeness, and since it is still 
possible that the light quark couplings are highly suppressed, making this the dominant 
production mode. The production/azimuthal angle distribution becomes 

72r . r ,»>2 

(18) 



oc sin 2 6 



d(4>e + (j)?) d cos 



(i) sin ( 2 cos (^ + 



5 One way to picture the construction of (4>i + <f>e) is as follows. Reflect the t in the tt production plane, 

thereby flipping the sign of its (f> angle. Then take the difference in (j> as before. 

6 This modulation is generically inducSdnbyTionzero rat effects (c/. Eq. 

53, 



7 Note that the Landau- Yang theorem 
gluon fusion. Implicitly, we only consider color octet resonances 

9 



54| forbids production of a color-singlet spin- 1 resonance through 



We see that the modulation persists, but it is now factorized from the production angle and 
is consequently twice as large when is integrated out. It also has the opposite sign. While 
this change in sign and magnitude could certainly lead to confusion over the true value 
of gL/gn, the production angle distribution is nonetheless still quite distinct, allowing an 
experimentally-driven determination of whether a Jbcam = or a Jb ca m = ±1 interpretation 
is appropriate. 



C. Spin-2 

Like spin-1, spin-2 particles can couple with different magnitudes to left and right-handed 
tops (and other fermions), though the leading interaction is now dimension-five, 

c oc (h{^v)i»t L ) + ^Lhr (t R (^uh,t R ) , (19) 

where h is a traceless and symmetric rank-2 tensor field, and A is the heavy scale that 
controls the strength of the operator. (This coupling is for a color-singlet field, but for 
color-octet we need only a trivial color matrix insertion.) Once again, we can build the 
traditional observables and see that they tell us nothing about the relative sign between left 
and right-handed couplings. (The double-cosine distributions and cosx distributions are 
identical to spin-1.) Just as with the other spins, the azimuthal angle distributions preserve 
this information. 

Quite analogous to spin-1, the top production and decay angle distributions are correlated 
with each other. To get the equivalent of Eq. [T6J, we must replace 

1 + cos 2 6 -> 1 - 3 cos 2 6 + 4 cos 4 e 

sin 2 6 -> - (l - 5 cos 2 6 + 4 cos 4 O) . (20) 

The coefficient of the modulating term now goes through a zero at cosG = 0.5. Integrating 
over production angles, the modulation is diluted with respect to spin-1 by a factor of three 
(and flipped in sign), 

« 1 + iQ 2 sin(2 f )cos(* + «. (2D 

While the production angle distribution is potentially a giveaway that the resonance is spin- 
2, the reduced (4>e + 4>i) modulation will make discrimination of the top couplings much more 
difficult. However, given adequate statistics, focusing on the region of central production 
would allow a much larger modulation to be seen. 

A spin-2 resonance can also be produced in gluon fusion, also through dimension-five 
operators. This can in principle lead to both Jbcam = ±2 and Jbcam = 0, but the latter 
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production mode is in fact absent at dimension-five I£| We therefore exclusively focus on 
■/beam = ±2. (We include Jbeam = in our discussions in Appendix (CJ as it may still arise 
from yet higher-dimension operators.) The production angle no longer factorizes out, and 
we are left with a joint production/azimuthal angle distribution much like Eq. [TBI but with 
an additional overall factor of sin 2 0. Since the production is biased more centrally, where 
the modulation is largest, the 0-integrated modulation is now enhanced, 

d l J r:t\ K 1 - I Q) 2 sm(2Ocos(0, + 0,). (22) 

Similar to the spin-1 case, gluon fusion leads to larger modulations with opposite sign. Again, 
distinguishing between qq and gg (as well as spin-1 in any production mode) is possible in 
principle by studying the distribution. 



III. OBSERVING AZIMUTHAL CORRELATIONS AT THE LHC 

We now turn to the question of whether the azimuthal modulations discussed in Section HT1 
are actually observable at the LHC, and in what cases we might be able to discriminate be- 
tween different resonances. We focus on the dileptonic mode, which has the largest azimuthal 
correlation effects amongst ti decay modes, but which also has the lowest branching frac- 
tion and, more pressingly, the fewest number of visible decay products. Given the presence 
of two experimentally-invisible neutrinos, it is fair to ask 1) whether the resonance can be 
identified, and 2) whether the azimuthal observables discussed above can be reliably recon- 
structed. We address these two concerns in turn, in the context of simple LHC simulations 
that include showering/hadronization, gaussian energy smearing, and basic analysis cuts. 
As described in more detail in the subsections below, we find that minimalistic strategies 
work quite well, and that the modulations are very robust. 



A. Simulation and basic event selection in the dileptonic mode 

We generate 6-body dileptonic resonance and background samples at leading-order, with 



full angular correlations, using MadGraph/MadEvent v4.4.51 [55j and its topBSM 23] 



For color singlet spin-2, the interaction h fJlV Tr(G fJ ' a G'^), where G is the gluon field-strength, contains no 
coupling to the Jbcam = state on-shell. We could instead try to couple through the parity-odd interaction 
h flv r Tr{G lia G'^) 1 but this vanishes identically when hu V is traceless. Considering an octet resonance, our 
only options are Tr(/i Mi , {G^ a , G v a }) or Tr(/i M i,[G MQ , G^]), each of which also lead to vanishing amplitude 

for Jbeam = 0. 
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modelJ§ We subsequently shower /hadronize with PYTHIA 6.4 56] and cluster the resulting 



particles into jets with Fast Jet v2.4.2 [57] without detector simulation. We then ap- 
ply gaussian energy smearings a(E)/E to the jets (5.6GeV/E © 1.25^GeV/E © 0.033), 
electrons (0.02), and muons (0.05 x y/E/TW). 

For very energetic tops, the lepton may be separated from the 6-jet by a fairly small 
angle, possibly failing ordinary isolation criteria and subjecting us to semileptonic heavy 
flavor backgrounds. Nonetheless, these leptons can probabl y b e reliably discriminated using 



tracker-level "mini- isolation" criteria with shrinking cones 58|, |59| , or through a number of 



other simple means [60j, |6l|. We accept into our analysis events with exactly two opposite 



signed leptons (electrons or muons) that pass mini-isolation as described in [58j], and with 
Pt > 20 GeV and \r)\ < 2.5. This set of cuts (or some modest modification thereof) should 
be adequate to reduce heavy flavor backgrounds to a very sub dominant level, while keeping 
the majority of the signal. 



We reconstruct jets using the anti-fcr algorithm [62| with R = 0.4, keeping those with 
Pt > 50 GeV and \r)\ < 2.5. No fe-tags are applied. Each lepton is then paired with the 
hardest jet which satisfies A(f)(j,£) < tt/2 and m(j,£) < 200 GeV@ Cases with too few or 
non-unique pairings are discarded. 

Our missing energy reconstruction is very naive, and hopefully conservative: we simply 
balance the summed of the two leptons and the two "6-jets." 

We find that reducible backgrounds from £ + £ _ +jets (including taus) and Vy + W^ - +jets 
are quite small throughout our analysis after all cuts, each at the 1-10% level with respect 
to the irreducible dileptonic tt0 In what follows, we explicitly keep only the irreducible 
background, which is by far our biggest concern. 



9 We do not include interference effects between the resonance and Standard Model in the ensuing analysis, 
since we expect them to be largely washed-out when integrating over the peak and taking the absolute 
value \4>i + 4>i\. Nonetheless, we have explicitly checked a small set of realistic, fully interfered samples 
with 15% resonance width, pure vector couplings to light quarks (to maximize the persistence of the 
interference given the chirality-averagcd initial state), and vector or axial couplings to tops. We find 
that with S/ B ~ 1 in the resonance region, our azimuthal distributions are consistent with an incoherent 

average of the signal and background. 

10 For genuine boosted top decays, the 6-jets are quite narrow, and their directions should be easy to 

accurately measure using the tracker. We therefore assign no error beyond energy smearing in the recon- 
struction of the jet-lepton invariant mass, nor in any of our subsequent kinematic reconstructions. 

11 To reduce the former background to this level, we require a simple Z-veto that rejects same-flavor lepton 

pairs with m(£ + ,£~) = [75, 105] GeV. This has a very minor effect on our signal, but order-of-magnitude 
effect on £ + ^~+jets. 
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B. Reconstruction of m t i 



The most efficient way to separate our signal from the continuum tt background is to 
measure m t i and isolate the resonance region. However, since we will investigate spin corre- 
lations in the dileptonic channel, we face an inevitable complication due to the two missing 
neutrinos. Of course, by the time a dileptonic analysis becomes sensitive to the resonance, 
we should already know the resonance's mass (and perhaps also its width) quite well from 
measurements in Z+jets and all-hadronic modes. So if we can reconstruct m t t, even approx- 
imately, then we already know where to look. 

The two neutrinos represent six lost kinematic degrees of freedom, only two of which are 
recovered through the measurement of $ T . Complete kinematic reconstruction is nonetheless 
possible by demanding on-shell tops and Ws on the two sides of the event, which fix the 
other remaining four degrees of freedom through four constraints. For energetic top pairs, 
the possible kinematic ambiguity regarding which 6-jet belongs to which lepton is not a 
major issue, and solving for the complete system essentially reduces to finding the zeros of a 
quartic polynomial [631. There are always four complex solutions, of which either zero, two, 
or four might be real] 12 ! 

This type of direct solution approach is not without its dangers. The continuum tt 
background at the LHC rises steeply toward small m^, and even a modest rate of misrecon- 
struction can "avalanche" low-mass backgrounds into our high-mass signal region. Indeed, 
the full kinematic solution can be highly sensitive to the precise values of visible kinematic 
variables, as well as to the detailed choice of top and W masses used to find the solutions!^! 
The situation is particularly difficult for energetic tops, as all of the decay products, including 
the neutrinos, approximately fall on a single axis in the transverse plane, and the neutrinos 
are roughly back-to-back. Accurate event reconstruction then demands exceptionally good 
measurement of missing energy components. 

The alternative is to consider reconstructions which are safer from the perspective of 
overestimating the mass of backgrounds, but which do so by actively attempting to un- 
derestimate the true mass. A well-known example (advocated for dileptonic tt resonance 
reconstruction in 



651 ] ) is the cluster transverse mass, 



Mki = (jp 2 TYiB + ™lis + #r) - (^ V is + It) 2 , (23) 

with the "vis" subscript referring to the 4-vector sum of the two leptons and two 6-jets. 
This represents the smallest m t i consistent with the measured ft T , without any on-shell 

12 A more realistic experimental analysis might instead solve the system using a multivariate likelihood 

approach (see, e.g., |64|). allowing all of the measured transverse obscrvables (]£ T , Pt{U)i Pr{ji)) to vary, 



approximately within their errors, along with the unknown neutrino 3-momcnta. This always gives real 

solutions by construction, and multiple solutions can be ranked or weighted by their likelihoods. 
13 Even the natural Brcit-Wigner scatter in the true event- by- event top and W masses produces sizable tails 

in the reconstructed m t j at parton-level. 
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criteria for the tops and H^s0 We can also try to improve upon this a bit by utilizing the 
fact that the neutrinos should be roughly collinear with the visible products when the tops 
are boosted. To this end, we construct a kind of minimal single neutrino, which under- 
reconstructs the complete vv subsystem, by setting Pt{v) = $t and setting t] to match 
that of the closest lepton in 0. (Setting rj to match the rapidity of the closest b + I system 
yields nearly equivalent results.) In addition, we could consider the M e g- variable used in 
supersymmetry searches, which simply adds up the scalar pxs of all visible products and 
]Et- We could also choose to be even more minimalistic, and just take the invariant mass 
of the visible component of the event @ 

To determine a reasonable method for our present work, amongst the many available 
options, we study a handful of kinematic reconstructions. In addition to the four under- 
complete reconstructions just described, we also consider two methods based the complete 
reconstruction using mass constraints, for which we must solve a quartic equation. While 
the quartic will practically never yield exactly one real solution, and often yields zero real 
solutions, we find that a simple and conservative workaround is to take the real components 
of the neutrino momenta returned by the solutions, and pick the configuration that yields 
the smallest m t i. We also include the strategy of Bai and Han 52] , which discards solutions 
that give the neutrinos too large a fraction of the scalar p? of the event, or which have 
imaginary components that are too large. Complex solutions are then made real by taking 
the complex norm of the reconstructed top and anti-top 4-momentum components, and all 
acceptable solutions are included with equal weight. (Complex solutions are only counted 
once.) Equivalently for our purposes, we just pick one solution at random. 

In Figs. [1] and [2j we show the results of our reconstructions on narrow 1 and 2 TeV spin- 
1 vector-coupled resonances and on the SM continuum at the 14 TeV LHC, including; an 
additional reference reconstruction where the neutrino 3- vectors are perfectly measured ] 16 " 17 ! 
Roughly speaking, there is a tradeoff between how tightly we attempt to sculpt the core of the 

14 More specifically, this is realized by a kinematic configuration where the neutrinos are exactly collinear in 

3D (the precise energy-sharing is irrelevant), are aligned with ]£ T in the transverse plane, and have rapidity 



matching that of the visible cluster. For generalizations of this method to cases where the missing 



particles 



GG 



67 1 . 



are massive (they should again be set 3D-collinear, but all moving with the same velocity), see 
15 Pe rhap s the most minimalistic strategy, simply applying a hard cut to the pt of the leptons, was used 



in 68|. This is fairly inefficient, so we do not pursue it here. However, in cases where the resonance 
has enough cross section, this may be a clean and effective way to get a reliable signal sample. We note 
that for highly boosted tops, the lepton px is largely uncorrelated with the azimuthal angles about the 
production axis, and we suspect that our correlation measurements will continue to work even with such 

a manifestly kinematically-biased cut. 

16 Axial resonances look nearly identical, while purely RH chiral resonances tend to have slightly harder 

distributions, and purely LH slightly softer. 

17 The reconstruction with perfectly measured neutrinos displays a pronounced low-mass tail, partially due 
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FIG. 1. Various approximations to m t i as applied to narrow 1 and 2 TeV vector resonances. 
Neutrino reconstructions include perfect (grey, scaled vertically by 1/2), Mt c \ (blue), minimal r\- 
collinear (red), visible-only (black), M e g (orange), minimal real-part quartic (green), and Bai-Han 
quartic (purple). 




FIG. 2. Various approximations to m t i as applied to continuum Standard Model tt at LHC14- 
Neutrino reconstructions include perfect (grey), Mnp c \ (blue), minimal rj-collinear (red), visible-only 
(black), M e ff (orange), minimal real-part quartic (green), and Bai-Han quartic (purple). 
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LHC7, M = 1 TeV 



LHC14, M = 2 TeV 




m 2.2 




FIG. 3. Significance improvement versus signal efficiency, using optimized mass windows and 
transversity cuts, for a 1 TeV narrow resonance at LHC7, and for a 2 TeV narrow resonance 
at LHC14- Neutrino reconstructions include perfect (grey, scaled to (l/2)es/ ' \[£b), Mt c \ (blue), 
minimal -q-collinear (red), visible-only (black), M e g (orange), minimal real-part quartic (green), 
and Bai-Han quartic (purple). (The 2 TeV curves were generated using a background sample with 
m t i > 800 GeV at parton-level, to improve statistics. Therefore the regions near 100% efficiency 
are not fully accurate. They would have a much steeper slope in reality, and the overall SIC 
levels on the rest of the plot would be higher. However, even with modest cuts, sensitivity to this 
generator-level cutoff largely disappears.) 



signal peak and how often we overshoot the true mass. To assess what this means in terms 
of discriminating signal from background, we run a simple cut optimization study, scanning 
over possible mass windows. Inspired in part by the observations of 52|, we also scan over a 
cut on the event's "transversity," defined as the ratio of either M^, M e Q, or min({p^(6, /)}) 
to the reconstructed mass. In particular, this type of cut aids in removing a substantial 
contribution from gluon fusion backgrounds, which tend to peak at low production angles. 
It also significantly cuts down the unphysical high-mass broadening of the background for 
the quartic reconstructions. 

Some representative results for the narrow resonance optimizations are presented in Fig. |3J 

to additional neutrinos lost in semileptonic i?-hadron decays, and partially due to FSR of gluons off of 

the tops. The shape can be improved with the inclusion of extra jets into the reconstruction, though we 

have erred on the conservative side and not done so. 
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LHC7 


Mass/Transversity Cuts 


-.selection 

£ s 


select ion-)~cuts 




1 TeV, Narrow 
1 TeV, 15%- Width 


M Tcl = [750,1025] GeV, M cG /M Tcl > 0.65 
M Td = [700,1450] GeV, M eS /M Td > 0.65 


0.38 
0.33 


0.24 
0.19 


21 fb 
24 fb 


2 TeV, Narrow 
2 TeV, 15%- Width 


M Td = [1600,2100] GeV, M cS /M Td > 0.50 
M Td = [1425,2925] GeV, M cS /M Td > 0.75 


0.58 
0.40 


0.20 
0.09 


0.07 fb 
0.17 fb 



TABLE I. Optimized mass window and transversity cuts, and the corresponding signal efficiencies 
and background cross sections, for the 7 TeV LHC. Signal efficiencies are for a vector-like spin-1 
resonance in the dileptonic mode, after basic selection and then after application of the cuts. 



LHC14 


Mass/Transversity Cuts 


^.selection 
£ S 


sclection+cuts 


OB 


1 TeV, Narrow 
1 TeV, 15%- Width 


M Tcl = [725,1000] GeV, M cS /M Tcl > 0.75 
M Tcl = [700,1300] GeV, M cS /M Tcl > 0.75 


0.35 
0.32 


0.19 
0.16 


179 fb 
253 fb 


2 TeV, Narrow 
2 TeV, 15%- Width 


M Tcl = [1475,2050] GeV, M eS /M Tcl > 0.60 
M Tcl = [1425,2425] GeV, M eB /M Tc \ > 0.60 


0.56 
0.48 


0.27 
0.19 


5.6 fb 
7.5 fb 



TABLE II. Optimized mass window and transversity cuts, and the corresponding signal efficiencies 
and background cross sections, for the 14 TeV LHC. Signal efficiencies are for a vector-like spin-1 
resonance in the dileptonic mode, after basic selection and then after application of the cuts. 



using the SIC ("significance improvement characteristic," i.e. £s/\Z^b) visualization of 48 ]. 
In practice, we find that M c g- usually works best for the transversity cut, and we use it for all 
of our displayed SIC curves except for that of M eS itself, for which we use mm({p T (b, /)}). 
We find that the consistently most powerful mass estimator is Mt c \, in combination with 
the cut on M eS /M Tc i. We therefore use these to define the signal regions for our subsequent 
azimuthal angle analysis, maximizing esjyfsB. Tables HI and ITTl contain the results of the 
optimization. 

Of course, any conclusions drawn from such a simple study can only be taken as sugges- 
tive, but it does indicate that attempts to make detailed reconstructions of the neutrinos 
might be counterproductive, unless cleverer techniques are brought to bear0 It is also not 
clear that the maximum ss/s/sb point is truly the best to maximize the significance of an 
angular correlation measurement. We make this choice merely for want of a better guide, 
and, in particular, lack of reliable information about possible systematic errors. 

18 We emphasize that the undcrpcrformance of the full quartic solutions with respect to Mt c \ in this study 

is not just an artifact of our somewhat pessimistic missing energy reconstruction. The conclusion remains 



unchanged if we switch to perfect f, T measurement and turn off energy smearings. 
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C. Measurement of azimuthal observables 



In the previous subsection, we investigated six different methods for reconstructing the 
complete dileptonic ti system. Here, we will now determine to what extent these methods 
can be used to extract azimuthal correlations. 

One issue which we must address at this point is how to define the production plane in 



the reconstructed tt CM frame, given that the beams are typically not exact 
in this frame. A common tactic is to follow the Collins-Soper construction 



y back-to-back 
, which defines 



an effective beam axis by taking the difference between the two unit vectors pointing along 
the beams. We take a slightly simpler route, which we find leads to some improvements 
in the amplitudes of our azimuthal modulations. We simply bring the ttsystem to rest in 
lab coordinates via a rotation-free boost, and then use the lab's native beam axis for the 
construction of the production plane. 

In Figs. H] through El we show the reconstructed distributions of {<t>t — <i>t) and \4>(, + 4>(\ for 
a variety of narrow gg-initiated spin-0 resonances and gg-initiated spin-1 resonances, and for 
the SM continuum at the 14 TeV LHC We apply our narrow resonance cuts from Table UTI to 
our resonance samples, and broad cuts to the SM samples (which slightly improves statistics). 
The reconstructions are the same as those discussed in the previous subsection, though we 
now omit M e s, which does not offer a fully 3D picture of the event. The figures exhibit all 
of the expected behavior discussed in Section |TTJ for all of the reconstruction techniques!^! 
Some are clearly better than others, but interestingly the differences are usually not very 
dramatic. The clear exception is the minimal real-part quartic solution, which tends to 
introduce spurious modulations in (<f)£ — (pi). However, the remaining reconstructions are 
all quite good. In particular, we preserve nearly the entire modulation using only visible 
particles (the two leptons aad two Hets)§ We take this as a strong hint that these angles 
may be very straightforward to reconstruct in reality, at least in the case of somewhat fast- 



19 We do acquire some small \<j>t + <\>i \ modulations for the scalars and (4>i — (f>i) modulations for the vectors, 
up to the 10% scale for 1 TeV resonances. These are mainly due to some biases incurred by our selection 
cuts and jct-lcpton pairing strategy, rather than genuine hclicity interference. Nonetheless, the effects of 

these biases also show some correlations with the resonance couplings. 

20 The visible-only 4>is are actually very closely related to angular variables which can be constructed in 

the top frame, and which are associated with a well-defined analyzing power of about 80%: instead of 
directly measuring the lepton direction, first boost the b + 1 composite system to rest (without rotation) , 
and then measure the lepton direction. In our visible construction, we instead immediately boost the 
b + / system into lab- frame without first bringing it to rest in top- frame, but the resulting 4>i is nearly the 
same. It is even possible to accurately approximate this using purely geometric constructions in the lab 
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1 TeV Spin-0~| 




2 TeV Spin-0~| 




FIG. 4. Distributions of ((pi — (pi) and \(pi + (pi\ for 1 and 2 TeV spin-0 resonances: pure scalar 
(solid), pseudoscalar (long- dashed), and mixed a = ir/4 (short-dashed). Neutrino reconstructions 
include perfect (grey), Mt c \ (blue), minimal rj-collinear (red), visible-only (black), minimal real-part 
quartic (green), and Bai-Han quartic (purple). 



frame. For example, each of the 6-jets can be taken to define an independent z-axis and production plane, 
with 4>£ 's measured with respect to these coordinate systems. We thank Stephen Parke for suggesting the 
possibility of using lab-frame angles. 
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1 TeV Spin-1 




2 TeV Spin-1~| 




FIG. 5. Distributions of (<j)£ — <pi) and \<p£ + (f>(\ for 1 and 2 TeV spin-1 resonances: pure vector 
(solid), axial-vector (long- dashed), and left-chirality (short- dashed). Neutrino reconstructions in- 
clude perfect (grey), Mrp c \ (blue), minimal r]-collinear (red), visible-only (black), minimal real-part 
quartic (green), and Bai-Han quartic (purple). 

moving tops where the visible products can be cleanly separated into hemispheres. The fact 
that the semileptonic top decay is roughly 2/3 visible is clearly a tremendous help here, as 
the vector-sum of its visible activity gives us an adequate approximation of the true top 
direction of flight. Moreover, our cuts bias us away from kinematic configurations where the 
neutrinos are dominating. 
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1 TeV Standard Model 




2 TeV Standard Model 




FIG. 6. Distributions of {(pi — (pi) and \<j>i + 4>e\ for Standard Model it at LHC14 in 1 and 2 
TeV (broad) mass windows. Neutrino reconstructions include perfect (grey), Mj- C i (blue), minimal 
n-collinear neutrino (red), visible-only (black), minimal real-part quartic (green), and Bai-Han 
quartic (purple). 

A detailed understanding of the observability of these modulations in the presence of 
background and systematic uncertainties is beyond the scope of this paper. However, we 
can make some rough estimates for the timely question of discriminating a true vector from 
an axial- vector. For this purpose, we can think of the \<f>i + <pe\ distribution with only 
two bins, i.e. as an asymmetry of roughly ±25%. From a purely statistical perspective, 
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with no background, these two cases could be discriminated at 3<r level if the error on the 
asymmetry can be brought down to the 15% level. This would require less than 100 events. 
Given that, after cuts, we accept only about 20% of the dileptonic signal, and hence 1% 
of the total it resonance signal, this kind of analysis becomes statistically possible when 
the LHC accumulates 0(10,000) it pairs from the resonance. We would therefore require 
cross sections at the pb-scale for the LHC7 (assuming 5-10 fb _1 ), and at the 100 fb-scale 
for the long-term LHC14. Interestingly, for £+jets, which has six times higher branching 
fraction but 40% as much modulation (if we correlate with the 6-quark on the hadronic side), 
the required cross sections would naively not be so different. On the other hand, different 
factors come into play, in particular the need to identify 6-(subJjets and to possibly run 



a hadronic top-tagger (as in, e.g., [60 . 



70M76(| . and reviewed in [77J), but also presumably 



larger efficiency and better resonance peak reconstruction. 

Of course, if S/B < 1, the discrimination becomes more challenging. Besides diluting the 
signal, the Standard Model continuum also exhibits its own vector-like (0£ + 0£) modulation. 
Accurate extraction of the signal modulation would therefore require a well-controlled back- 
ground subtraction. Assuming this can be done, then models with S/B ~ 1 could probably 
be discriminated at 3a with at least a few hundred dileptonic signal and background events 
after cuts. From Table [J we see that the SM backgrounds to a 1 (2) TeV analysis after 
10 fb" 1 at the LHC7 would be roughly 200(1) events. We are background-limited only for 
resonances at or slightly above 1 TeV, requiring multi-pb cross sections. Heavier resonances 
are rate limited. For 100 fb" 1 at LHC14, even 2 TeV resonances will be background-limited, 
requiring low pb-scale for the signal. While we have not investigated even higher resonance 
masses, these again quickly become rate-limited at the 14 TeV machine, with 100 fb yielding 
adequate statistics to discriminate vector from axial-vector in the dileptonic mode. 

Given these estimates, what might we expect to learn about axigluon-like models appli- 
cable to the Tevatron top Apb anomaly? The size of the asymmetry (central value near 
50% for nan > 450 GeV [l|) and the lack of major deviations from the SM prediction for 



da/drritt |78l . 1 791 ] favor a heavy resonance with large couplings. At the same time, LHC 



searches for dijet resonances and contact interactions |80l . |8l| constrain the couplings to 
quarks. While we do not undertake a full scan of the parameter space, we can use the 
analysis of [lo| to identify a reasonable example of a phenomenological model, consistent 
at 68% CL with the Tevatron anomaly and having modest impact on dcr/dm t f. a 1.5 TeV 
resonance with purely axial couplings to light quarks and top quarks of g\ = (0.7) g s and 
g l A = (— 3.0)g s o Assuming that b L also couples with the same strength as the tops, the 
width is about 18%, and the branching fraction to tops is close to 57% (to light quarks it is 
about 14%). This resonance easily evades the exclusion limits on universal axigluons and on 

21 The contribution to Afb hi the m t j > 450 GeV bin would be about 27% at parton- level. Combined with 

the NLO contribution from the Standard Model, this brings the asymmetry to within la of the measured 

value. 
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more general dijet resonances in the recent ATLAS search 80j, is safe from quark contact 
interaction constraints^ and is below the current CMS limit on top resonances 82]. At 
LHC7, a x BR(ti) is about 4 pb, and the optimal Mt c \ window should be highly signal 
dominated. A 5-10 fb _1 data set should therefore be enough to establish a pure axial in- 
terpretation over a pure vector interpretation to at least 3a for this model. Clearly, then, 
azimuthal decay correlations could potentially have a nontrivial part to play in disentangling 
the full story behind the Tevatron Afb anomaly, even for the early phase LHC. 

It is also worth noting that our analysis could still apply to cases without a well-defined 
resonance peak. For axigluon-like resonances well above 1 TeV, matching the size of the 
Tevatron anomaly requires large couplings, and therefore typically quite large width (T/M > 
20%). This may simply show up as a broad excess over the Standard Model at high mass. 
However, neglecting for the moment interference with the Standard Model, none of our 
results were actually sensitive to where we looked on the resonance's Breit-Wigner curve. 
Indeed, while a broad resonance signal may be difficult to isolate on-peak, it can come to 
dominate the high-mass tail of ti production. With a full 14 TeV LHC at 0(100) fb -1 
luminosity, this tail will be very well-populated by the Standard Model alone. For example, 
the region above 2 TeV should contain 0(10,000) central QCD-induced events, and this 
number may be substantially enhanced by the resonance. 

Still, the final fate of such an analysis becomes more model-dependent, as interference 
with the Standard Model could be non-negligible. This affects not only the total rate, but 
also the coefficient of cos(<^ + (pi), in proportion to the product of the light quark and top 
quark vector charges. Notably, this interference is washed-out if either the top quarks or the 
light quarks are mostly axially coupled. See Appendix |B] for a more detailed discussion. 



IV. CONCLUSIONS 

The discovery of a resonance in the ti invariant mass spectrum at the LHC would represent 
the beginning of a major advance in our understanding of TeV-scale physics. In the event 
of such a discovery, determining the precise nature of this resonance will become one of the 
highest priorities of the LHC experimental program. Fortunately, the very distinctive decays 
of top quarks will give us the opportunity to address rather detailed questions about the 
resonance's couplings. Unlike the case of a resonance decaying directly into stable particles 
such as electrons, different ti helicity channels coherently interfere as they decay into the 
same set of 6-body final-states. The final decay distributions, correlated between top and 
antitop, therefore carry information about signs/phases in the tops' chiral production matrix 
elements. 

22 We have explicitly checked that the contribution to central dijet production in each of ATLAS's mass 
bins, including full interference and t-channcl exchanges, is below the stated uncertainties. 
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We have argued in this paper that the tops' azimuthal decay distributions about their 
production axis provide a sensitive and experimentally robust probe of these correlations. 
As is usual for it spin correlations, the charged lepton (or down-type quark) exhibits the 
largest effects. Spin-0 resonances lead to 60% modulations in (<fre — <^), with a phase offset 
that directly represents the resonance's CP phase. Spin-1 and spin-2 resonances lead to 
modulations in (<p£ + <^), with a signed magnitude that depends on the ratio of chiral 
couplings, (30%) x 2 / {ql/qr + Qr/ql) when integrated over production angles (two times 
larger when restricting to highly central production). While the spin-0 CP phase has long 
been a topic of investigation, the spin-1 (2) chiral coupling ratio is largely missed in common 
top spin correlation variables such as the double-cosine distribution and the 3D lepton 
opening angle. We have found that single-differential azimuthal sum/difference distributions 
contain all of the relevant information contained in these more common variables, and more. 

We have also found that these azimuthal distributions are straightforward to reconstruct, 
without detailed information on neutrino kinematics. We established this point in the dilep- 
tonic mode, in which variables sensitive to spin correlations usually require complicated 
reconstruction techniques. In contrast, azimuthal distributions can be well-approximated 
using only visible kinematics. The major nontrivial task then becomes isolating the reso- 
nance region by reconstructing Our own (admittedly limited) investigation suggests 
that simple reconstructions such as Mt c \ may be superior to more complete reconstructions 
that utilize on-shellness criteria on the Ws and tops. However, this conclusion could still 
change in a more detailed experimental study, or with better reconstruction techniques. 

Since the modulations are large in the dileptonic mode, measuring them should not require 
a very big event sample. For example, 0(100) events should be adequate to discriminate 
a pure vector resonance from an axial-vector resonance at better than 3a level, assuming 
S/B > 1. For favorable choices of parameters, a color-octet axial- vector resonance could 
explain the Tevatron Aps anomaly, be discovered at the early LHC, and have its mostly- 
axial coupling to tops confirmed not long afterwards. This would be a dramatic interplay 
between the capabilities of the two colliders. 

While we have focused on dileptonic it, the option of using the £+jets mode will also 
be important to explore in more detail. Extracting the smaller modulations may require 
greater control over systematics, but the much higher statistics and much better kinematic 
reconstruction afforded by this channel could easily offset this. In particular, for resonances 
well above 2 TeV, £+]ets may be the only viable option, even for the long-term LHC, due 
to the smaller cross sections involved. 

We have also limited ourselves to simulations based on leading-order it production and 
decay, with additional radiation and kinematic rearrangements provided by the PYTHIA 
virtuality-ordered parton shower. Understanding the genuine NLO behavior of the azimuthal 
correlations, processed through realistic event reconstructions, remains as an interesting open 



task which can potentially be addressed by modern calculations |83l-l85|. 
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Finally, while the focus of this work has been characterization of physics beyond the 
Standard Model, we note that azimuthal decay angles could also be useful for detecting the 
spin correlations of boosted or semi-boosted tops from pure QCD (cf. the modulations in 
Fig. [H]), and serve as a particularly clean indication that the top quark is spin- 1/2. 
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Appendix A: Basic Definitions and Spin-0 

In Section [III we gave various formulas for dT over different portions of the angular phase 
space, including the production angle of the ti pair in the ti CM frame and the decay angles 
for the two leptons. In these appendices, we present these formulas in a more complete form, 
displaying all angular correlations. We also go beyond the chiral m t = limit. Though this is 
likely to be a reasonable approximation for resonances at 1 TeV or higher, finite-mass effects 
are not always totally negligible. We parametrize them with factors of r = 2mt/M (the 
inverse Lorentz boost of the tops) and \J\ — r 2 (the tops' velocities). However, we continue 
to ignore velocity-suppression factors that arise purely from ti phase space, which in any 
case can be factorized out of the total resonance decay rate. We also neglect interference 
terms with the SM. These effects tend to be small on-peak. There are, of course, cases where 
the interference is nonetheless more important than the finite-mass corrections. We outline 
out the effects of interference in more detail below. 

Recalling Section [Til we define our coordinate system as follows. We first find the ti 
production axis in the ti CM frame. The angle that this makes with the beam axis is called 
G, and we average over configurations 6 f>7r-6. We then perform rotation- free boosts to 
bring both tops to rest. We define a common z-axis using the original t direction. To define 
the y-axis, we take either of the two unit vectors perpendicular to ti production plane, and 
then define the x-axis to make a right-handed 3D coordinate system. 

Before proceeding, we note that all of the equations that we give can apply to arbitrary 
pairs of top decay products. For example, in the £+jets channel we might correlate the 
lepton with the 6-jet from the hadronic side. Starting with the full angular distributions for 
two leptons, upon changing to the angular variables for the particles of interest, we need 
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only multiply each sin and cos of the new polar angle by the appropriate spin- analyzing 

power: 1 for leptons, -0.3 for neutrinos, -0.4 (+0.4) for 6-quarks (W-bosons), and 0.5 for the 

softer of the two W decay products in the top rest frame r 3 H 24 l We subsequently drop the 'f 

subscript from all of the angular variables, as an implicit reminder that we are no longer 

restricted to dileptonic mode. 

The angular formulas for spin-0 decay were given in Section III Al to zeroth order in r. 

Restoring the full kinematic dependence, we get 

d A T -/ 

- oc 1 + cos 9 cos 9 — sin 9 sin 9 ( cos(2a) cos (0 — 0) 
ail ail V 

— y/l — r 2 sin(2o;) sin(0 — 

i 

■r 2 (l + cos(2a)) (l + cos 9 cos 9 — sin 9 sin 9 cos(0 — 



dr 



2 

oc ( 1 — r 2 cos 2 a 



d(0- 

— ^— J ^ [ cos(2a) — r 2 cos 2 a\ cos(0 — 0) — y/l — r 2 sin(2a) sin(0 — 0) J ,(A1) 

where 2a is the relative phase between the ii^R and i^ti, couplings. Values of a besides and 
7r/2 signify CP violation. The upper formula is the "master" leptonic angular distribution 
from which one can derive the dependence on whatever angular variable one wishes. The 
lower equation describes the azimuthal correlations, which we found useful for identifying a 
heavy particle's relative couplings to different top chiralities regardless of its spin. In this 
case, the dependence is in (0 — 0). For spin-0, (0 + 0) is a physically meaningless quantity, 
though we will see below that spins 1 and 2 have their leading modulations in it. Formulas 
for the traditional distributions in cos 9 cos 9 and cos x are given in Section [TTJ 

In addition to the finite-boost corrections, we can also consider the effects of interference 
with the Standard Model. First, we note that these effects vanish for a scalar produced in 
qq annihilation, since to very good approximation the QCD contribution only takes place 
between quarks of opposite helicity (same spin). For gg fusion, the interference vanishes 
at high boost, as the QCD contribution dominantly proceeds from same-spin gluons to 
same-spin (chiral) tops. The interference is also suppressed relative to the "resonance- 
squared" contribution since the interference term crosses a zero on-peak. Integrating the 
interference contribution across the peak (accounting for imperfect cancellation due to falling 
PDFs and/or asymmetric signal window), the total relative suppression goes as r(r/M) 2 . 
In the case of a color-singlet scalar with generic couplings to G a ^ u G a ^ u and G a ^ v G a ^ u , the 

23 Note that care should be taken if, for example, we decide instead to define the z-axis by following the 

direction of the semileptonic top in £+jets. In that case, the formulas for the i? + +jets events stay the 

same, but for the £~+jcts events they require replacements 9 — > tt — 9 and 9 — > n — 9. 

The analyzing powers for particles other than the lepton are only formally correct under the assumption 

that the phase space of the remaining particles is fully integrated. In principle, they could change after 

accounting for detector acceptance and analysis cuts. 

26 



interference in the (<ft — <p) modulation is controlled by the latter coupling for production 
nearly perpendicular to the beamline (and swamped by the quasi-singular pure QCD rate 
for production near the beamline). While the resonance-squared contribution exhibits a 
modulation cos(0 — <ft + 2a), the interference comes in with sin(0 — <p + a), for example 
signalling CP-violation when the resonance couples to tops as a pure scalar (a = 0) and 
to gluons as a pure pseudoscalar. For a 1 TeV resonance of 20% width, we estimate that 
the size of the interference contribution relative to the resonance- squared is less than 20% 
if S/B > 1/4. This relative contribution decreases with increasing mass, decreasing T/M, 
or increasing; s/bE 



Appendix B: Spin-1 

We give the tops' coupling to a spin-1 particle in Eq. [TTJ Since the interactions involve 
same-chirality fields, the left- and right-handed couplings are independent in the absence of 
parity. We parametrize the different chiral couplings as follows, 

g L = g cos f 

g R = gsm£. (Bl) 

Just as with the scalar case, including m t ^ in the resonance decay will modify the 
differential decay rate. Since the equivalence of chirality and helicity breaks down with 
finite fermion mass, all possible spin pairings contribute, though only pj" and \\ remain in 
the limit r — > 0. 

As mentioned above, producing a spin-1 resonance with zero angular momentum along 
the beam axis requires gluon fusion, which only occurs via a higher-dimension operator. We 
therefore begin with formulas for the likely more dominant scenario of Ji, cam — il via qq 
annihilation. 

Including the dependence on the production angle and all decay angles (choosing one 
particle each from the t side and t side), we get 
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Quite generally, we can infer that the relative contribution of interference within a mass window of width 



O(r) centered on the resonance is at most (r/M)/y S/B, with a typically 0{\) coefficient that depends 
on the specific process. In our case, this coefficient carries an additional suppression factor of r. 
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+ sin 2 sin 9 sin 9 (cos cos + sin(2£ ) sin sin 



(B2) 



From this distribution, we can integrate and/or change variables to get a function of whatever 
angular variables we want. We start by demonstrating the claim made in Section [III that 
for a sufficiently light resonance, one can determine the extent to which it is vector or axial 
with the traditional cos 6 cos 9 distribution. We can see the difference by setting £ = 7r/4 



or 37r/4 and integrating the other angles in Eq. 
prefactor, we get 



Dividing out by a common numerical 
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Firstly, we note that this reproduces the values listed in Table II of Ref. 23[. We see 
that the ratio of modulating to constant term for the vector only approaches that of an 
axial resonance in the chiral limit. This difference at nonzero r allows one to separate the 
two cases, though for M above 1 TeV, this requires fitting the 2D distribution to percent- 
scale precision. Lastly, we note that the decay amplitude vanishes for an axial-vector at it 
threshold. This is simple to understand with parity. An axial vector is even under P, while 
the exchange of ti in the final state is odd. For the fermion pair at rest, there is no orbital 
contribution to the wavefunction that can make up the difference, so the amplitude shuts 
off. We give the distributions in (0 + 0), (0 — 0) in Eq JB5l 

For completeness, we include the analogue of Eq. IB2I for Jb cam = 0, as would arise for a 
color-octet produced in gluon fusion. Production of such a state is usually suppressed due to 
the need for higher-dimension interactions, but could dominate if dimension-four couplings 
to light quarks are small. The 5D distribution in leptonic and production angles is 

d 5 r Jh = 

° beam 
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From Eqs. IB2I and |B4 | we derive azimuthal distributions for polarizations ±1 and 0, 

(l + ^r 2 (3sin(20 - l)) - Cp (£)* (sin(2£) + ^r 2 (l - sin(2£)) ] cos(0 + 



dTj, =p 

'-'beam - 1 

d((f) + 0) 
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d(0-0) 



dVj — =P oc (l + ^r 2 (3sin(20 - l)) - \r 2 (£)* (l + sin(20) cos(0 - 0), (B5) 



where (p = (1/2, —1) for P = (±1, 0), respectively. It is straightforward to see why the 
different polarizations have the same (0 — 0) dependence. The (0 — 0) modulation arises 
from mass-suppressed "wrong-helicity" (tl or |f) tt production. Once 6 is integrated out, 
this spinless configuration of tops cannot display any residual dependence on the resonance 
polarization. 

Our discussion thus far has focused on the resonance-squared contribution to the total 
rate, but there is also interference with SM tt production. This is a small effect on-peak 
for narrow resonances, and even for broad resonances if the couplings are large compared 
to QCD. Still, in a full analysis, one should account for such effects, especially as they 
may be numerically important relative to the subleading r 2 corrections in e.g. Eq. IB5I Let 
us consider in detail the case of qq annihilation into a spin-1 resonance. The interference 
comes out proportional to the light quarks' vector charges (having averaged over forward and 
backward directions as well as incoming quark chiralities) , and consists of two terms. One 
is a cos(0 + 0) modulation, which simply tracks the total cross section interference with the 
SM, and is proportional to the top's vector coupling. This passes through a zero on-peak, as 
the resonance propagator goes 90-degrees out of phase with the gluon propagator. Thus, we 
get a sizable interference effect on-peak only if the light quark and top quark charges each 
have relatively large vector components, and if the off-peak contributions in our resonance 
mass window are highly imbalanced. The second interference term is proportional to the 
top axial charge, and introduces a sin(0 + 0) contribution to the modulation. This term 
is suppressed by T/M, and would have been averaged out in our analysis of Section [III C\ 
which only looked at |0 + 0|. (Note that a sinusoidal modulation component in (0 + 0) does 
not signal CP- violation, as it would in (0 — 0).) 



Appendix C: Spin-2 

The spin-2 case shares many features with spin-1, including separate couplings to left 
and right chirality fields. We give the lagrangian for this case in Eq. [19j Once again, we 
parametrize the relative strengths of the chiral couplings with an angle £ (Eq. IB 1 [) . As 
mentioned in Section [TT], qq and gg production give rise to Jbcam = ±1 and Jb cam = ±2 spin 
states, respectively, and each proceed via dimension-five operators. We therefore should not 
prejudice one over the other in general, though gg production tends to benefit from higher 



parton luminosities in the case of lighter spin-2 resonances (less than about 2 TeV 23| ) , and 
qq for heavier resonances. The 5D angular distributions are 
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Just as with spin-1, the cos 9 cos 9 distribution lets us distinguish in principle the difference 
between vector-like and axial-like couplings. For all spin-2 polarizations the relevant formulas 
are 
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d cos 9 d cos 9 

where V and A are respectively the cases of left-right symmetric and antisymmetric cou- 
plings. We see just a slight modification in the vector term compared to spin-1, but once 
again having M > 1 TeV requires us to fit the distribution to within a few percent to 
discriminate between the two cases. 

To avoid needing such precision and a full 2D parameter fit, we can integrate the lepton 
polar and production angles to find the distributions in 0,0. The leading term in r only 
depends on + 0. We get 
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where (p = (—2/3, 1/6, 1) for P = (±2, ±1, 0), respectively. Once again, we see that the 
(0 — 0)- dependence is independent of resonance polarization (cf. discussion under Eq. IB5j) . 

The production of a spin-2 particle with Jb ca m = requires an operator with dimension 
greater than five, i.e. beyond that of the other polarizations. It is for this reason that we 
did not list any equations for this state in Section [Til However, we give here the differential 
decay rate with respect to leptonic and production angles, from which one can obtain the 
0-polarization portion of Eq. f 
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